Interatomic collisions in a tightly confined Bose gas 
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We discuss pair interatomic collisions in a Bose gas tightly confined in one (axial) direction and 
identify two regimes of scattering. In the quasi2D regime, where the confinement frequency coo 
greatly exceeds the gas temperature T, the scattering rates exhibit 2D features of the particle motion. 
At temperatures T ~ huio one has a confinement-dominated 3D regime, where the confinement can 
change the momentum dependence of the scattering amplitudes. We describe the collision-induced 
energy exchange between the axial and radial degrees of freedom and analyze recent experiments on 
thermalization and spin relaxation rates in a tightly (axially) confined gas of Cs atoms. 
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I. INTRODUCTION 

Collisional properties of ultra-cold gases strongly con- 
fined in one direction attract a great deal of interest since 
the start of active studies of spin-polarized atomic hydro- 
gen. In the latter case the interest was related to recom- 
bination and spin relaxation collisions and to elastic scat- 
tering in the (quasi) 2D gas of atomic hydrogen adsorbed 
on liquid He surface (see [jy for review). The discov- 
ery of Bose-Einstein condensation in trapped alkali-atom 
clouds stimulated a progress in evaporative and op- 
tical cooling and in trapping of neutral atoms. Present 
facilities make it possible to (tightly) confine the motion 
of particles in one direction to zero point oscillations. 
Then, kinematically the gas is 2D, and the only differ- 
ence from the purely 2D case is related to the value of the 
interparticle interaction which now depends on the tight 
confinement. Thus, one now has many more opportuni- 
ties to create (quasi)2D gases. In the recent experiments 
with optically trapped Cs J5|-|| about 90% of atoms are 
accumulated in the ground state of the harmonic oscilla- 
tor potential in the direction of the tight confinement. 

In this paper we consider a Bose gas tightly confined 
in one (axial) direction and discuss how the axial con- 
finement manifests itself in pair elastic and inelastic col- 
lisions. We identify two regimes of scattering. At tem- 
peratures T -C huo (luq is the axial frequency) only the 
ground state of the axial harmonic oscillator is occupied, 
and one has a quasi2D regime. In this case, the 2D char- 
acter of the relative motion of particles at large sepa- 
ration between them, manifests itself in a logarithmic 
energy dependence of the scattering amplitude. For a 
negative 3D scattering length a, we observe resonances 
in the dependence of the elastic scattering rate on a. This 
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is quite different from the 3D case where the scattering 
rate always increases with a 2 . The presence of these res- 
onances in quasi2D follows from the analysis given in || 
and finds its origin in increasing role of the 2D kinematics 
of the particle motion with increasing ratio |a|/Zo, where 
Iq = (h/muio) 1 / 2 is the axial extension of the atom wave- 
function, and m the atom mass. 

At temperatures T ~ htoo we have a confinement- 
dominated 3D regime of scattering, where the 2D charac- 
ter of the particle motion is no longer pronounced in the 
scattering process, but the axial confinement can strongly 
influence the energy (temperature) dependence of the 
scattering rate. Treating collisions as three-dimensional, 
the wavevector p of the relative motion of colliding atoms 
does not decrease with T. The atoms undergo zero-point 
oscillations in the axial direction and this corresponds to 
p ~ 1/Zo- If the 3D scattering amplitude is momentum- 
dependent at these p, which is the case for \a\ > Z , then 
the temperature dependence of the elastic collisional rate 
becomes much weaker. This means that for a large 3D 
scattering length the tight axial confinement suppresses a 
resonant enhancement of the collisional rate at low ener- 
gies. In many of the current experiments with ultra-cold 
gases one tunes a to large positive or negative values by 
varying the magnetic field and achieving Feshbach res- 
onances |]ro|-]l4|. In the unitarity limit (\a\ — ► oo) the 
3D elastic cross-section is a — 8ir/p 2 and the rate of 
3D elastic collisions strongly increases with decreasing 
temperature. The tight confinement of the axial motion 
makes the scattering rate practically temperature inde- 
pendent at T ~ hujQ. We obtain a similar suppression 
of resonances for inelastic collisions, where the resonant 
temperature dependence in 3D is related to the energy 
dependence of the initial wavefunction of colliding atoms. 
We analyze the Stanford and ENS experiments on elastic 
HH and spin relaxation || collisions in a tightly axially 
confined gas of cesium atoms and discuss the origin of sig- 
nificant deviations of the observed collisional rates from 
the 3D behavior. 

We develop a theory to describe the collision-induced 
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energy exchange between axial and radial degrees of free- 
dom of the particle motion. We establish selection rules 
for transitions between particle states in the axial har- 
monic potential and calculate the corresponding transi- 
tion amplitudes. This allows us to consider tempera- 
tures T > hioo and analyze thermalization rates in non- 
equilibrium clouds. In the Stanford and ENS experi- 
ments these clouds were created by means of degenerate 
Raman sideband cooling which effectively leads to 
a gas with different axial (T z ) and radial (T p ) tempera- 
tures. After the cooling is switched off, the temperatures 
T z and T p start to approach each other, and ultimately 
the gas reaches the equilibrium temperature. At suffi- 
ciently low T only a few axial states are occupied and 
the temperature dependence of the corresponding ther- 
malization rates should deviate from the 3D behavior, 
thus exhibiting the influence of the axial confinement on 
the scattering process. We calculate the thermalization 
rates and establish the conditions under which this influ- 
ence is pronounced. 

The minimum energy exchange between the radial and 
axial degrees of freedom of two colliding atoms is equal 
to 2Tiu!q. This follows from the symmetry of the inter- 
atomic potential with respect to simultaneous inversion 
of the axial coordinates of the two atoms, which ensures 
the conservation of parity of their wavefunction under 
this operation. Accordingly, the sum of two (axial) vi- 
brational quantum numbers can be changed only by an 
even value. The rate of energy transfer from the ra- 
dial to axial motion is proportional to the difference be- 
tween the radial and axial temperatures AT = T p — T x , 
if they are close to each other. As the total energy of 
colliding particles should exceed 2haj in order to en- 
able the energy transfer, the rate of this process at 
temperatures T < Tiloq becomes exponentially small: 
AE cx ATexp(— 2fkJo/T). Due to the presence of the 
energy gap Hloq in the excitation spectrum of the axial 
harmonic oscillator, the heat capacity of the axial degree 
of freedom is dE z /dT z ~ exp(— huio/T), which leads to a 
thermalization rate AT/ AT cx exp(— TujJq/T). This ex- 
ponential temperature dependence shows that the ther- 
malization is suppressed at very low temperatures. One 
can deeply cool the axial motion, but radially the cloud 
remains "hot" on a very long time scale. 



II. 2D SCATTERING PROBLEM 



The quantity /(q, <fi) is the scattering amplitude, q is the 
relative momentum of the atoms, and <fr the scattering 
angle. Note that f(q,<fi) in Eq.([j]) differs by a factor of 
— y/8nq from the 2D scattering amplitude defined in [fL5| . 

Similarly to the 3D case, the scattering amplitude is 
governed by the contribution of the s-wave scattering if 
the relative momentum q satisfies the inequality qR e <C 
1, where R e is the characteristic radius of interaction. In 
the case of alkali atoms, the radius R e is determined by 
the Van der Waals tail of the potential U(p) and ranges 
from 20 A for Li to 100 A for Cs. The s-wave scattering 
amplitude is independent of the scattering angle <fi. The 
probability a(q) for a scattered particle to pass through a 
circle of radius p per unit time is equal to the intensity of 
the scattered wave multiplied by 2irpv, where v — 2hq/m 
is the relative velocity of colliding atoms. From Eq.(|l]) 
we have 



a(q) = —\f(q)\ 2 . 
2 m 



(2) 



The velocity v is equal to the current density in the inci- 
dent wave of Eq. ([!]). The ratio of a(q) to this quantity is 
the 2D cross section which has the dimension of length: 



a(q) = \f(q)\ 2 /4q. 



(3) 



For the case of identical bosons Eqs. (||) and (||) have an 
extra factor 2 in the rhs. 

The quantity a(q) is nothing else than the rate con- 
stant of elastic collisions at a given q. The average of 
a(q) over the momentum distribution of atoms, multi- 
plied by the number of pairs of atoms in a unit volume, 
gives the number of scattering events in this volume per 
unit time. 

For finding the s-wave scattering amplitude one has 
to solve the Schrodinger equation for the s-wave of the 
relative motion of colliding atoms at energy e — h 2 q 2 /m: 



h 2 



u(p) 



h 2 q 2 

tps{q, p) = 4>s(q, P)- (4) 



At distances p ^> R e the relative motion is free and one 
can omit the interaction between atoms. Then the solu- 
tion of Eq . (|4|) , which for qp I gives the partial s-wave 
of ip(p) P); takes the form 



ip s (q,p) = Jo(qp) 



H {qp), p<i? e 



(5) 



First, we discuss the purely 2D elastic scattering in 
pair collisions of ultra-cold atoms interacting via a short- 
range potential U(p). At interparticle distances p — > oo 
the wavefunction of colliding atoms is represented as a 
superposition of the incident plane wave and scattered 
circular wave 



f{q, 4>)\ 



a l <lP 



8irqp 



(f) 



where Jo and Hq are the Bessel and Hankel functions. 

On the other hand, at distances p< 1/<j one can omit 
the relative energy of particles in Eq. (Q) . The resulting 
(zero energy) solution depends on the momentum q only 
through a normalization coefficient. In the interval of 
distances where R e <C p <C l/q, the motion is free and 
this solution becomes ip s cx \n(p/d), where d > is a 
characteristic length that depends on a detailed shape of 
the potential U (p) and has to be found from the exact so- 
lution of Eq.M) with q = 0. This logarithmic expression 
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serves as a boundary condition for ip s (q, p) (|J) at qp <C 1 , 
which immediately leads to the scattering amplitude fl5} | 



2n 



\n{l /qd*)+m/T 



(6) 



where d* = (d/2) expC and C « 0.577 is the Euler con- 
stant. 

It is important to mention that the condition qR e <C 1 
is sufficient for the validity of Eq. (||) . This equation also 
holds for the case of resonance scattering, where the po- 
tential U(p) supports a real (or virtual) weakly bound 
s-level. In this case the spatial shape of ip s (q, p) at dis- 
tances where R e -C p -C l/q, is the same as the shape 
of the wavefunction of the weakly bound state. This 
gives <i» = h/^/meo, where £o is the binding energy. 
We thus have the inequality <i* 3> R e , and the quan- 
tity qd* in Eq.(|6|) can be both small and large. The 
rate constant aia) peaks at q = 1/d, and decreases as 
1/[1 + (4/7r 2 )ln [qd*)\ with increasing or decreasing q. 
Note that the 2D resonance is actually a resonance in 
the logarithmic scale of energies. The decrease of a by 
factor 2 from its maximum value requires a change of 
energy e = h 2 q 2 /m by factor 20. 

For qd* <C 1 one may omit the imaginary part in 
Eq.(^|), and the scattering amplitude becomes real and 
positive |Tq] . The positive sign of f(q) has a crucial 
consequence for the mean-field interparticle interaction 
in purely 2D Bose gases. In the ultra-cold limit where 
qR e <C 1, the scattering amplitude is related to the en- 
ergy of interaction in a pair of particles (coupling con- 
stant g). For a short-range potential U(p), the energy 
of the mean-field interaction in a weakly interacting gas 
is the sum of all pair interactions. In a uniform Bosc- 
condensed gas the coupling constant g for condensate 
atoms is equal to the amplitude of scattering (with an 
extra factor ft 2 /m for our definition of /) at the energy 
of the relative motion e — h 2 q 2 /m — 2p, where p is the 
chemical potential (L8|. Hence, we have 



.9 = — J(<lc 
m 



2irh 2 



1 



m ln(l/<7 c d*) 



> 0; q c d, < 1, (7) 



where q c — ^/2mp/ft is the inverse healing length. In a 
dilute thermal 2D gas, due to the logarithmic dependence 
of / on q, the thermal average of the mean-field interac- 
tion leads to the coupling constant g = (ft 2 / m) f (qx) , 
where qT = V mT /ft is the thermal momentum of parti- 
cles. At sufficiently low temperatures, where qxd* <C 1, 
we again have g > 0. 

Thus, in an ultra-cold purely 2D gas the coupling con- 
stant for the mean-field interaction is always positive in 
the dilute limit and, hence, the interaction is repulsive. 
This striking difference from the 3D case is a consequence 
of the 2D kinematics. For low energies, at interparticle 
distances /) » i? e , the (free) relative motion of a pair 
of atoms is governed by the wavefunction ip g oc ln(p/d). 
The probability density IV'sl 2 of finding two atoms at a 



given separation increases with p as the condition p > d 
is always reached, unless the atoms have a bound state 
with energy e — > [d — > oo). This means that it is favor- 
able for particles to be at larger p, i.e. they repel each 
other. 



III. SCATTERING IN AXIALLY CONFINED 
GEOMETRIES. GENERAL APPROACH 

In this Section we discuss elastic scattering of atoms 
(tightly) confined in the axial (z) direction, assuming 
that the motion in two other (x, y) directions is free. We 
analyze how the scattering is influenced by the confine- 
ment and calculate a complete set of scattering ampli- 
tudes corresponding to collision-induced transitions be- 
tween particle states in the confining potential. We still 
call this scattering elastic as the internal states of atoms 
are not changing. 

For a harmonic axial confinement, the motion of two 
atoms interacting with each other via the potential V(r) 
can be still separated into their relative and center-of- 
mass motion. The latter drops out of the scattering prob- 
lem. The relative motion is governed by the potential 
V(r), together with the potential Vh(z) — uo 2 z 2 /4 origi- 
nating from the axial confinement with frequency luq. For 
the incident wave characterized by the wavevector q of 
the motion in the x, y plane and by the quantum number 
v of the state in the potential Vh(z), the wavefunction 
of the relative motion satisfies the Schrodinger equation 



A + V(r) + V H {z)- — - 

m 2 



<A(r) = e^(r), (8) 



where e = H q /m + vTwJq. 

The scattering depends crucially on the relation be- 
tween the radius of interatomic interaction R e and the 
characteristic de Broglie wavelength of particles A e . The 
latter is introduced qualitatively, as the motion along 
the z axis is tightly confined. Accounting for the zero 
point axial oscillations one can write A e ~ ftj\J me, with 
i = e+huj n /2. We will consider the ultra-cold limit where 



Ae > Re 



(9) 



Eq.(^|) immediately leads to the inequality qR e <C 1, as 
the de Broglie wavelength for the motion in the x, y plane 
is ~ l/q. For small v the harmonic oscillator length 
Iq = (ft/muso) 1 / 2 plays the role of the axial de Broglie 
wavelength of atoms. Therefore, the ultra-cold limit (j|) 
also requires the condition Iq 3> R e - For large v, the 
axial de Broglie wavelength is ~ Iq/ \Jv and, according to 
Eq. (^|) , this quantity should be much larger than R e . 

Under the condition qR e <C f, the scattering ampli- 
tudes are determined by the contribution of the s-wave 
for the motion in the x, y plane. In the case of iden- 
tical bosons, the s-wave scattering requires even values 
of v and v' as the wavefunction tp should conserve its 
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sign under the transformation z — » — z. The quantum 
numbers v and v' should be even also for distinguishable 
particles. Otherwise at distances of interatomic interac- 
tion, r < R e , the wavefunction ip will be small at least 
as R e /lo, ensuring the presence of this small parameter 
in the expressions for the scattering amplitudes. 

The scattering amplitudes corresponding to transitions 
from the initial state v (of the relative motion in the 
potential Vh(z)) to final states v' are defined through 
the asymptotic form of the wavefunction ip at an infinite 
separation p in the x, y plane: 



ip(T)m<p v (z)e 



fvv<(s)ip v <{z)i 



(10) 



where ip v (z) and <p v '{z) are the (real) eigenfunctions of 
the states v and v' . For each of the scattered circular 
waves the value of the momentum q v i follows from the 
energy conservation law fi. 2 g?//m = e — hio^v' > 0. 

Relying on the condition (Q), we develop a method that 
allows us to express the scattering amplitudes through 
the 3D scattering length. At interparticle distances 
r 3> R e the relative motion in the x, y plane is free, 
and the motion along the z axis is governed only by the 
harmonic oscillator potential Vh{z). Then, the solution 
of Eq.(|8|) with V(r) — can be expressed through the 
Green function G e (r,r') of this equation. Retaining only 
the s-wave for the motion in the x, y plane, we have 

TP(r) = <p v (z)J (qp)+A l/ G e {r,0), (11) 

and the expression for the Green function G e (r, 0) reads 



G E (r,0)=$^v(2)M0)x- 



iH$-\q v ,p)/±; ql,>0 
K (\q v ,\p) /2tt; q 2 v , <0 



(12) 



Here the summation is also performed over closed scat- 
tering channels for which q 2 , < 0. The function Kq(x) = 
(iTt/2)Ho(ix) and it decays as y/n/2x exp (—x) at x ^S> 1. 
Thus, for p — > oo the terms corresponding to the closed 
channels vanish. Then, comparing Eq.(nTI) at p — > oo 
with Eq. (|l0|) , we find a relation between the scattering 
amplitudes and the coefficients A v : 



f vvl = -A v (p v >(0)9(e - hu! v'), 



(13) 



where 9 is the step function. 

The condition Iq ^> R e ensures that the relative motion 
of atoms in the region of interatomic interaction is not 
influenced by the axial (tight) confinement. Therefore, 
the wavefunction ^(r) in the interval of distances where 
R e <C r <C A e , differs only by a normalization coefficient 
from the 3D wavefunction of free motion at zero energy. 
Writing this coefficient as (p v (0)r), we have 



ip(r) w tp u (0)rj(l - a/r). 



(14) 



Eq.(pT[) serves as a boundary condition for ip(r) ( pi] ) at 
r -> 0. 



For r — > 0, a straightforward calculation of the sum in 
Eq.@ yields 



G s (r,0) 



1 



1 



( £ 



2(2tt) 3 / 2 Z \2hLu 
where the complex function w(x) is given by 



w(x)- 



- lim 



N N 
In — 

7T 



(2j-l)!! 



N 

(2/)!! 



bx(j— x— iO) 



(15) 



(16) 



With the Green function (|15|), the wavefunction (|llj) at 
r — > should coincide with i/j(r) (|T4|). This gives the 
coefficient 

V = TT-T-, ; ; (17) 



1 + (a/V^lo)w(e/2huJo) 

and provides us with the values of the coefficients A v . 
Then, using Eq.(|l^) and explicit expressions = 
[1 / 2-kIq) 1 / \v — l)\\/y/iA, we immediately obtain the scat- 
tering amplitude foo(s) and express all other scattering 
amplitudes through this quantity: 



/oo(e) = 47r^o(0)a?7 



2V2tt 



lo/a + (l/V2ir)w(£/2huio) 



(18) 



f vv ,(e) = P vv ,f m {e)6{e - hu; v)9(e - fo^i/), (19) 



where 



= ^(ow(o) = { v -iy.\{ v '-iy} (2Q) 



^(o) 



One can see from Eqs. (yj|) and (|l9| ) that for any transi- 
tion v — ► v 1 the scattering amplitude is a universal func- 
tion of the parameters cl/Iq and e/Tuoq. The quantity 
P uv i in Eq. ( p^ ) is nothing else than the relative probabil- 
ity amplitude of having an axial interparticle separation 
\z\ <C Iq (in particular, \z\ < R e ) for both incoming [y) 
and outgoing (z/) channels of the scattering process. It 
is thus sufficient to study only the behavior of foo(e). 

We emphasize the presence of two distinct regimes of 
scattering. The first one, which we call quasi2D, re- 
quires relative energies e -C Tiujq. In this case, the rel- 
ative motion of particles is confined to zero point os- 
cillations in the axial direction, and the 2D kinematics 
of the relative motion at interatomic distances p > Iq 
should manifest itself in the dependence of the scatter- 
ing amplitude on ej2hiOQ and cl/Iq. In the other regime, 
at energies already comparable with Tiujq, the 2D kine- 
matics is no longer pronounced in the scattering pro- 
cess. Nevertheless, the latter is still influenced by the 
(tight) axial confinement. Qualitatively, the scattering 
amplitudes become three-dimensional, with a momen- 
tum ~ l/lo related to the quantum character of the axial 
motion. Thus, we can say that this is a confinement- 
dominated 3D regime of scattering. With increasing the 
relative energy to e > fiuo, the momentum is increas- 
ing to y/ms/h and the confinement-dominated 3D regime 
continuously transforms to ordinary 3D scattering. 
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IV. QUASI2D REGIME 

In the quasi2D regime, due to the condition e <C Tiojq, 
the incident and scattered waves have quantum numbers 
v = v = for the motion in the axial harmonic po- 
tential Vh{z). The relative energy e = h 2 q 2 /m and the 
inequality qlo <C 1 is satisfied. In this case Eq.dlrf) gives 



w(e/2fiujo) = In (BTwJq / ire) + in, 



(21) 



where B w 0.915. Then our equatio n (|lq ) recovers 
Eq.(ll) of ref. 1 9], obtained in this limit []20|. 

Using Eq.(pi) we can represent /00(e) (|lq) in the 2D 
form (^j), with 



d* = (d/2)expC = ■s/tt/B l exp (- yV^ h/a). (22) 

This fact has a physical explanation. Relying on the 
same arguments as in the purely 2D case, one finds that 
in the interval of distances where Iq <C p <C l/q, the 
wavefunction tj) oc tpo{z) hi (p/d). On the other hand, for 
p > (o we have ^(r) = (po(z)ip s {p), where ip s is given 
by the 2D expression (pi ) with f(q) = /00(e)- This fol- 
lows from Eqs. (|Tl|) -(|l5| ), as all closed scattering chan- 
nels iy' 7^ 0) in Eq. (|12|) for the Green function G e (r,0) 
have momenta \q v i | > 1/Iq and will be exponentially sup- 
pressed at p ^> Iq. Matching the two expressions for the 
wavefunction ip one immediately obtains the 2D equation 
(^). However, the parameter d* ( ^2| ) can be found only 
from the solution of the quasi2D scattering problem. 

We thus conclude that the scattering problem in the 
quasi2D regime is equivalent to the scattering in an effec- 
tive purely 2D potential which leads to the same value of 
d*. For positive a <ti Iq, this potential can be viewed as a 
(low) barrier, with a height Vo ~ % 2 a/mlo and radius Iq. 
Hence, in the case of positive a we have a small (positive) 
scattering amplitude, in accordance with Eqs. (^|) and 
(p2|). For a negative a satisfying the condition \a\ <C Iq, 
the effective potential is a shallow well which has a depth 
I Vo| and radius Iq. This shallow well supports a weakly 
bound state with an exponentially small binding energy 
£0, which leads to an exponentially large as follows 
from Eq. (p2|) . As a result, we have a resonance energy 
dependence of the scattering amplitude /oo at a fixed ra- 
tio a/lo, and a resonance behavior of /oo as a function of 
a/lo at a fixed e/hujQ. 

The resonance in the energy dependence of /oo is quite 
similar to the logarithmic-scale resonance in the purely 
2D case, discussed in Section [n]. The quasi2D resonance 
is also described by Eq.(||), where the length d* is now 
given by Eq.(^). As expected, the dependence of /oo on 
e is smooth. 

The resonance in the dependence of the quasi2D scat- 
tering amplitude on a/lo has been found and discussed 
in |5). Relying on the above introduced effective 2D po- 
tential for the quasi2D scattering, we can now explain 
this resonance on the same grounds as the resonance in 
the energy dependence of /oo- We will do this in terms 



of the relative energy e and the binding energy in the 
effective potential, £0 = fi 2 /md 2 oc exp(Zo/|a|). For 
e/eq — (qd*) 2 ^> 1, the scattering amplitude in Eq.(||) 
is real and negative. It increases in magnitude with de- 
creasing ratio e/£o, that is with decreasing q or Iq. In 
the opposite limit, where s/sq <C 1, the scattering ampli- 
tude is real and positive and it increases with the ratio 
£/eo- The region of energies £/eo ~ 1 corresponds to 
the resonance, where both the real and imaginary parts 
of /oo are important. The real part reaches its maxi- 
mum at £/£o = exp (— 7r), drops to zero at £/eo = 1, and 
acquires the maximum negative value for s/sq = exp7r. 
The dependence of Im/oo on e/eq is the same as that of 
the quantity |/oo| 2 - Both of them peak at e/eq = 1 and 
decrease with increasing or decreasing e/eq. 

Qualitatively, the picture remains the same for \a\ ~ ^o- 
In Fig.l we present the dependence of |/oo| 2 on E/2hojQ 
at cl/Iq equal to —1, 1, and 00. In the two last cases 
we always have e/eq <§; 1, and |/oo| 2 increases with £ at 
£ Uloq. For a/l = —1 we have the above described 
logarithmic-scale resonance in the behavior of |/oo| 2 - 



15- 



10- 




10"' 
e/2ftco 

FIG. 1. The function |/oo| 2 versus energy for a/lo — — 1 
(solid curve), a/lo — 00 (dashed curve) and a/lo = 1 (dotted 
curve) . 

The quasi2D resonance is much more pronounced in 
the dependence of the scattering amplitude on the pa- 
rameter cl/Iq. The reason is that /oo logarithmically de- 
pends on the particle energy, whereas the dependence on 
a/lo is a power law. For £ <g; Tiloq, Eq.(18) yields 



in 



16tt 2 



mi 



(V2nl Q /a + In {BTluq/tie)) 2 + tt 2 



(23) 



The quantity |/oo| 2 differs only by a factor of h/m from 
the rate constant of elastic collisions (see Eq.(pl)), and 
one can think of observing the resonance dependence of 
|/oo| 2 on a/lo m a n experiment. For example, one can 
keep £ (temperature) and cl>o constant and vary a by using 
Feshbach resonances. The resonance is achieved at a = 
— Iq In (BTuoa / ~ke) . This is a striking difference from the 
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purely 3D case, where the cross section and rate constant 
of elastic collisions monotonously increase with a 2 . 

In Fig. 2 we present |/oo| 2 versus a/lg at a fixed e/huiQ. 
In order to extend the results to the region of energies 
where the validity of the quasi2D approach is question- 
able, the quantity |/oo| 2 was calculated by using Eq.jlgl) 
for the scattering amplitude. The resonance is still visible 
at e/huJo = 0.06 and it disappears for e/Tiluq = 0.2. 

The obtained results allow us to conclude that for \a\ > 
lo the approximate border line between the quasi2D and 
confinement-dominated 3D regimes is e ~ e* = O.ITiljq. 
For | a | <C lo, as we will see below, the confinement- 
dominated regime is practically absent. 




FIG. 2. The function |/oo| 2 versus |a|/?o at various energies 
for a < (a) and a > (b). 

The output of kinetic studies in thermal gases is usu- 
ally related to the mean collisional frequency (the rate of 
interatomic collisions) f2 = cm, where a is the mean rate 
constant of elastic collisions, and n the gas density. In 
the quasi2D regime, the rate constant a follows directly 
from Eq.(|^), with twice as large rhs for identical bosons: 



broadens the resonances in Fig.l and Fig. 2. Due to the 
logarithmic dependence of /oo on the relative energy, the 
thermal average is obtained with a good accuracy if one 
simply replaces e by the gas temperature T. Thus, in 
order to observe the manifestation of the 2D features of 
the particle motion in their collisional rates one has to 
achieve very low temperatures T < O.ltkoo- 



V. CONFINEMENT-DOMINATED 3D REGIME 

In the confinement-dominated 3D regime, where e ~ 
Hujq, the axial confinement influences the scattering pro- 
cess through the confined character of the axial motion. 
In order to analyze this influence, we first examine the 
function w(e/2huj ) which determines the energy depen- 
dence of the scattering amplitudes. The imaginary part 
of w(x), following from Eq. (|l6|) , is equal to 



Inru;(x) = 7r 

3=0 



(2J-1)!! 
(2j)» 



2^ 



r(N + 3/2) 



(25) 



where [x] is the integer part of x. The function Imra(i) 
has a step- wise behavior as shown in Fig. 3. It is constant 
at non-integer x and undergoes a jump at each integer 
x, taking a larger value for larger x. With increasing 
x, the jumps become smaller and for i > 1 we have 
lmw(x) ps 2^/tFx. The real part of w(x) was calculated 
numerically from Eq.(|l6|) and is also given in Fig. 3. At 
any x we have |Reu>(a;)| < 1, except for narrow intervals 
in the vicinity of integer x. In each of these intervals 
the function Kew(x) logarithmically goes to infinity as x 
approaches the corresponding integer value. This is con- 
sistent with the step- wise behavior of lmw(x) : As one can 
see directly from Eq.(|l^), for x approaching an integer j 
the analytical complex function w oc In (J — x — iO). 
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where the symbol ( ) stands for the thermal average. Our 
numerical calculations show that the average of |/oo| 2 
over the Boltzmann distribution of particles only slightly 



FIG. 3. The functions Kew(x) (solid curve) and IniK;(x) 
(dashed lines). The dotted curve shows the function 
corresponding to the asymptotic behavior of Imw at large x. 
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The described behavior of the function w(e/2hu>o) has 
a direct influence on the scattering amplitudes. For 
e/2hu!o close to an integer j, the amplitude is small and 
it is equal to zero for e = 2hu>oj. This phenomenon 
originates from the fact that for e close to 2hu>oj, a 
new scattering channel opens (really or virtually). For 
this channel the momentum \qu'\ — \frn\e — 2hu>oj\/h is 
very small. Hence, at distances p <C |l/gy| the wave- 



function tp (11) will be determined by the contribution 
of this low-momentum term if p 3> Iq. This is clearly 
seen from Eqs. ( |l2| ) and ( |TT| ) and makes the situation 
somewhat similar to that in the quasi2D regime of scat- 
tering. In the latter case, the wavefunction ip ( |lT| ) at 
distances p <C 1/g is also determined by the contribu- 
tion of the low- momentum channel as long as p_3> ^o- 
Then, as follows from the analysis in Section IV, this 
wavefunction and the scattering amplitude f 00 behave as 
1/ In (Tlloq/e) in the limit e — > 0. In the present case, the 
wavefunction ip and the scattering amplitudes are small 
as l/ln(?iwo/|e — 2Tiu>oj\) for e — > 2frui j. 
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FIG. 4. The function |/oo| 2 versus energy. In (a) the pa- 
rameter a/lo = —1. In (b) the dashed curve corresponds 
to the unitarity limit (a/lo = oo), and the dotted curve to 
a/lo = 1- 

The energy dependence of |/oo| 2 on s/2hujQ for a/lo 
equal to —1, 1, and oo is displayed in Fig. 4. Outside nar- 
row energy intervals in the vicinity of integer s/2fiujo, the 
quantity |/oo| 2 is a smooth function of e. One can also 
see a sort of a step- wise decrease of |/oo| 2 with increasing 
e, originating from the step- wise increase of the function 
lmw(e/2hui ). For e/2huj close to integer values j > 



we find a fine structure similar in nature to the behavior 
of |/oo| 2 at e <§; hojQ. For a/lo = 1 and a/lo = oo there are 
narrow dips corresponding to the logarithmic decrease of 
|/oo| 2 as e — > 2huJoj, and for a/lo = — 1 these dips are 
accompanied by resonances. Note that the thermal dis- 
tribution of particles averages out this fine structure, and 
the latter will not be pronounced in kinetic properties. 

The difference between the confinement-dominated 3D 
regime and the ordinary 3D regime of scattering will 
manifest itself in the rate of elastic collisions (mean col- 
lisional frequency fi). For the Boltzmann distribution 
of particles, one can find this quantity by turning to 
the thermal distribution for the relative motion of collid- 
ing partners. Collision-induced transitions between the 
states of the relative motion in the axial potential Vh(z) 
are described by the rate constants 



ow(e) = (n/m)\f uv <(e)f 



(26) 



where the scattering amplitudes /„„/ are given by 
Eqs. (|l8| ) and (|l9|), and an extra factor 2 for identical 
bosons is taken into account. The collisional frequency 
51 = an, where n is the (2D) density, and the mean rate 
constant of elastic collisions, a, is obtained by averaging 
OLuv' ( |26] ) over the thermal distribution of relative ener- 
gies e and by making the summation over all possible 
scattering channels. We thus have 



£/w"-< £ »-M-f> (27) 



Here At = (2-kTi 2 /mT) 1 / 2 is the thermal de Broglie wave- 



length, e — h q /m + Tiujov, and the quantum numbers v 
and v' take only even values. The distribution function 
over v and q is normalized to unity, and the normaliza- 
tion coefficient A is equal to 



A = 2(l-cxp(-^o/? 1 )). 



(28) 



Note that the number of collisions per unit time and unit 
surface area in the x, y plane is equal to an 2 /2 = fin/2. 

The manifestation of the tight axial confinement of 
the particle motion in collisional rates depends on the 
relation between the scattering length a and the charac- 
teristic de Broglie wavelength A e ~ h/y/m(e + Tiljo/2) 
accounting for the zero point axial oscillations. For 
the scattering length satisfying the condition \a\ <C A E , 
the scattering amplitudes are energy independent at any 
e, except for extremely small energies in the quasi2D 
regime. This follows directly from Eqs. (pT[)-(^9[). The 
condition \a\ <C A e automatically leads to the inequal- 
ities \a\ <C Iq and \a\ <C fi/^/me. Hence, the func- 
tion w(e/2TiLJo) is much smaller than /o/|a|, unless e < 
hup exp (— /p/|a|) (see Eq.([2l|) and Fig. 3). Accordingly, 
Eq.@ gives J) = 1 and Eqs. (p|), @ lead to the scat- 
tering amplitudes 

f vv > = 4iratpv(0)<p u >{0)6(e ~ haj iy)9(e - hoj u'). (29) 
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The amplitudes ( |29| ) are nothing else than the 3D scat- 
tering amplitude averaged over the axial distribution of 
particles in the incoming (y) and outgoing (z/) scattering 
channels. From Eqs. ( |26| ) and (|2^) one obtains the same 
rate of transitions v — > v' as in the case of 3D scattering 
of particles harmonically confined in the axial direction 
and interacting with each other via the potential V(r). 
This is what one should expect, since under the condition 
\a\ <C A e the amplitude of 3D scattering is momentum 
independent. The integration over d 2 q in Eq.(p7|) leads 
to the mean collisional frequency 



fi= ^(47ra) 2 A^^ (0)(p 2 ;(0)exp ^_^O max{ ^ u , } 



Let us first show how the 3D result follows from our 
analysis at T > Filuq, irrespective of the relation between 
a and A E . At these temperatures the main contribution 
to the sum in Eq.(|27|) comes from e Tiloq and large 
v and v'. Accordingly, we can replace the summation 
over v and v' by integration. At energies much larger 
than huio the quantity y/v/l^ is nothing else than the 
axial momentum k z and we have e = 
these energies the function w(e/2Tiujo) in Eq.(|l8| 
its asymptotic form w ~ i \J 2tte / Tiujq . 
this immediately allows us to write 




takes 
Using Eq.@, 



\fvv' I — P vv 



l 2 (l+p 2 a 2 ) p ^ 



&3D m 
I 2 ' 



V < 



hujo ' 



Thus, in the case where \a\ <€. At, the tight confine- 
ment in the axial direction can manifest itself in the col- 
lisional rates only through the axial distribution of par- 
ticles and the discrete structure of quantum levels in the 
axial confining potential. The expression for the colli- 
sional frequency f2 can be reduced to the form 



(30) 



where the coefficient £ ranges from 1 at T <C hu>o to 2/w 
for T 3> hujQ. The condition \a\ <C A T is equivalent to 
\a\ <C Iq and qr\o\ "C 1, where qr — v 7 mT/h is the ther- 
mal momentum of particles. For T 3> hujQ, Eq.(|30|) gives 
the collisional frequency which coincides with the three- 
dimensional result averaged over the classical Boltzmann 
profile of the 3D density in the axial direction, ub{z): 



^3D — 



»(*) 



dz. 



(31) 



Here a is the 3D elastic cross section, and v is the rel- 
ative velocity of colliding particles. In other words, the 
quantity (l/2)an 2 = (l/2)f2n coincides with the number 
of 3D collisions per unit time and unit surface area in the 
x,y plane, given by (l/2)(a 3D v) J n 2 B (z)dz. 

From Eq.(|30|) we conclude that for \a\ <C lo the 
confinement-dominated 3D regime of scattering is not 
pronounced. At temperatures T < Kloq the collisional 
rate only slightly deviates from the ordinary 3D be- 
havior. This has a simple physical explanation. For 
\a\ <C Ay, treating collisions as three-dimensional we 
have fl ~ % r Ka 2 vn^£ ) . At low temperatures T < Tiloo the 
velocity v ~ %/mlo and the 3D density is n 3 o ~ ti/Iq. 
For T ^> hujQ we have v ~ {T/m) 1 / 2 and n 3 o ~ 
n(mcjQ/T) 1 / 2 . In both cases the "flux" vn 3 D ~ u>on, 
and there is only a small numerical difference between 
the low-T and high-T collisional frequencies. 

The ultra-cold limit @ assumes that the characteristic 
radius of interatomic interaction R e <^ Iq. Therefore, the 
condition |a| <gC Zo is always satisfied, unless the scatter- 
ing length is anomalously large (\a\ 3> R e ). Below we will 
focus our attention on this case, assuming that \a\ > No- 



where p = yj ' q 2 + k 2 is the 3D momentum of the relative 
motion, and a 3 £, = 87ra 2 /(l +p 2 a 2 ) is the cross section 
for the 3D elastic scattering. For large v and v' , Eq.(|20|) 



gives P vv > = (4/tt 2 



, and the integration over v' 



in Eq.(p7|) multiplies a 3 ^ by the relative 3D velocity v. 
Then, turning from the integration over v to the integra- 
tion over the axial momentum, we reduce Eq.(p7|) to 



nA T d 3 p 
(2tt) 3 



(a 3D v)Aexp 



h 2 p 2 
mT 



(32) 



and one can easily check that Eq.([52j) coincides with the 
three-dimensional result Q 3 d (|l]) . 

In the limiting case, where the thermal momentum of 
particles satisfies the inequality qr\ci\ 3> 1, we obtain 



167m fhuj \ , 

&3D = ; <fr ft > 1. 

m \ T J 



(33) 



In the opposite limit, where <7r|ft| -C 1, at temperatures 
T hujQ we automatically have \a\ -C Ay and, accord- 
ingly, recover Eq.(|30|) with ^ = 2/tt: 



16hn 



3D 



q T \a\ < 1. 



(34) 



As mentioned in Section for \a\ > Iq the approxi- 
mate border line between the quasi2D and confinement- 
dominated 3D regimes is e* rs OATlluo- In the temper- 
ature interval e* < T < Hcjq, the leading scattering 
channel will be the same as in the quasi2D case, that 
is v = v' = 0. However, the expression for the scatter- 
ing amplitude /oo is different. From Fig. 3 and Eq.(|l8|) 
one concludes that the real part of the function w can be 
neglected and the scattering amplitude takes the form 



/oo — 



2V2^ 



lo/a + i^7r/2 

Then, retaining only the scattering channel v = v 1 = 
Eqs. @ and (§?]) yield 

- exp (—huJo/T) 
Jo J 1 + ira 2 /2l 2 



(35) 
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The difference of Eq.(p5[) from the quasi2D result of 
Eqs. (|23l ) and ( |2~i| ) is related to the absence of the loga- 
rithmic term in the denominator. This follows from the 
fact that now we omitted the real part of the function w, 
which is logarithmically large in the quasi2D reg ime. 

It is worth to note that for lo ^> Jal, Eq.(p5|) is only 
slightly different from the 3D result (|34|). This is consis- 
tent with the above given analysis leading to Eq.([30|). 

On the other hand, for large \a\/lp the difference be- 
tween Eq.(|35|) and the 3D result (33) is significant. This 
originates from the fact that for a large scattering length 
a the 3D amplitude of scattering in the ultra-cold limit 
depends on the particle momenta. For a tight axial con- 
finement, treating collisions as three-dimensional, the rel- 
ative momentum of colliding particles at temperatures 
T < h^o is ~ 1/fo an( l it n0 longer depends on tem- 
perature. Hence, the scattering rate is quite different 
from that in 3D. Given these arguments, one expects a 
strongly pronounced confinement-dominated 3D regime 
of scattering if the ratio |a|/Zo ^S> 1. 
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FIG. 5. The quantity mQ/hn versus temperature. In (a) 
the parameter a/lo = — 1 (dashed curve), and a/lo = 1 (dot- 
ted curve). In (b) a/lo = 00 (unitarity limit). The solid 
curves in (a) and (b) show the 3D result (|3l]). The arrow in 
(b) indicates the lowest ratio T/%u>o in the Stanford and ENS 
experiments. 



This is confirmed by our numerical calculations for the 
temperature dependence of Vt from Eq.(p7|). In Fig. 5 we 
present the results for a/lo equal to — 1, 1, and 00. The 
largest deviation from the 3D regime is observed in the 
unitarity limit (a — > 00). From Fig. 5 we see that in 
the Stanford || and ENS || experiments performed in 
this limit one should have significant deviations of 
collisional rates from the ordinary 3D behavior. 



VI. THERMALIZATION RATES 

We will now discuss the collision-induced energy ex- 
change between axial and radial degrees of freedom of the 
particle motion in an ultra-cold Bose gas tightly confined 
in the axial direction of a pancake-shaped trap. It is as- 
sumed that the radial confinement is shallow and it does 
not influence the scattering amplitudes. In this geometry, 
using degenerate Raman sideband cooling, the Stanford 
HU and ENS |7j0] groups created Cs gas clouds with 
different axial (Tz) and radial (T p ) temperatures. Af- 
ter switching off the cooling, interatomic collisions lead 
to energy exchange between the axial and radial particle 
motion and the temperatures T z and T p start to approach 
each other. Ultimately, the gas reaches a new equilibrium 
state, with a temperature in between the initial T z and 
T p . The corresponding (thermalization) rate has been 
measured at Stanford || and ENS § and it provides us 
with the information on the regimes of interatomic colli- 
sions in the gas. 

The radial motion of particles is classical. Therefore, 
we will calculate the rate of energy exchange between the 
radial and axial degrees of freedom for a given value of 
the radial coordinate p and then average the result over 
the Boltzmann density profile in the radial direction. The 
latter is given by 



n(p) = n(0) exp 



? ? 
mui • p 

2T 



(36) 



where n(0) = moj 2 N / 2nT is the 2D density for p = 0, to 
the radial frequency, and N the total number of particles. 
Collision-induced transitions v — > v 1 change the energy of 
the axial motion by Tkooiy 1 — v). We will assume that in 
the course of evolution the axial and radial distribution 
of particles remain Boltzmann, with instantaneous values 
of T z and T p . Then the rate of energy transfer from the 
radial to axial motion can be written on the same grounds 
as Eq. (|27|) and reads 



- 1 iW'(e) I A exp 
m V ml p 1 , 



l/V ' 

h 2 q 2 hcjQis 



A 2 T d 2 q 



%ujq[v — v) x 



(37) 



where e — h 2 q 2 jm + Tiojqv, and the normalization coeffi- 
cient A depends now on both T z and T p . 







The radial energy of the gas is E p = 2NT p , and the 
axial energy is given by E z — Nhuj [exp (huj /T) — 1] . 
The time derivatives of these energies take the form 



E z 



Nfi 2 uj 2 f z 



E n = 2NT n 



(38) 



4T|sinh 2 (^ /2T z )' 

Given the initial values of T z and T p , Eqs.(|37]) and j3^ ) 
provide us with the necessary information on the evolu- 
tion of T z (t) and T x {t). 

For a small difference ST = T p — T z , these equations 
can be linearized with respect to ST. As the total energy 
is conserved, Eqs.(^) reduce to 



• 4/1, 4T 2 sinh 2 (^o/2T)\ 

ST= w{-2 + J 



(39) 



In Eq.(p7|) we represent the exponent as — [e/T + 
SThajQv/T 2 ) and turn from the integration over dq to 
integration over de. The zero order term of the expan- 
sion in powers of ST vanishes. The (leading) linear term, 
being substituted into Eq.([39|), leads to the differential 
equation for ST(t): 

St = -n th {T)ST, (40) 

where the thermahzation rate £lth(T) is given by 



th 



16nh 



n(0)A 2 T A (h 2 u 2 , A ^ 2 



2T 2 



4sinli (haj /2T) x 



poo 

J2(v-v') 2 cfe|jw( £ )| 2 ex P 



(41) 



The degeneracy parameter is n(0)A^ = N(hu>/T) 2 and 
it is small as the gas obeys the Boltzmann statistics. The 
normalization coefficient A is again given by Eq. (p8[) . 

The quantum numbers v and v' take only even val- 
ues and, hence, in order to change the state of the ax- 
ial motion one should have a relative energy e > 2Tujq. 
Therefore, at temperatures lower than huo the rate of 
transitions changing the axial and radial energy is cx 
exp (— 2TujJq/T). On the other hand, the axial energy 
E z cx exp (— Hloq/T) and thus the thermahzation rate 
Qth oc exp (— Tiujq/T). This can be easily found from 
Eq.(|l]) and shows that the quantum character of the 
axially confined particle motion exponentially suppresses 
the thermahzation process at temperatures T <C %u)q. In 
particular, this is the case for the quasi2D regime. 

In the most interesting part of the confinement- 
dominated 3D regime, where e* < T < %ujq, the energy 
exchange between the axial and radial motion of parti- 
cles is mostly related to transitions between the states 
with v' = and v = 2. The relative energy e should be 
larger than 2Tiujq and, at the same time, this energy is 
well below 4hu>o- Hence, the scattering amplitude /20(e) 
is determined by Eqs. ( |l8| an d (jl9| ) in which one can put 
w(e/2Hlu ) w i37r/2 (see Eq.(|25D and Fig. 3). This gives 



h 




and from Eq.([H]) we obtain 

n IQhwQ _ exp(— hajn/T) , , _ 

°» = 1 + (1 ~ ^JT)f- (42) 

The characteristic frequency f^o is given by 

flo = lu 2 N/uj . (43) 

At temperatures T ^> fnug, Eq.([ll]) leads to the 3D 
result for the thermahzation rate: 



3D 
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15tt 
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Q ; q T \a\ > 1, 



(44) 



Comparing Eqs. (|44| ) and ( (45|) with Eq. (f42|) one sees that 
VL t h should acquire its maximum value at T ~ hu>o- For 
\a\ ^ this maximum value is on the order of £lo/2ir. 

As one expects from the discussion in Section M, the 
difference of the thermahzation rate from Q^P is pro- 
nounced for large values of a. For example, in the 
unitarity limit Eq.(E3) gives n^P cx 1/T 2 , whereas 
in the confinement-dominated regime we have flth oc 
(l/T)exp(-huo/T). 

It should be emphasized that for any T, ujq, and a the 
ratio Qth/Qo depends only on the parameters T/Hloq and 
a/l^. This can be found directly from Eq.([i"l|). In Fig.6 
we present the temperature dependence of Qth in the 
unitarity limit, obtained numerically from Eq.( |4l|), and 
compare our results with the data of the Stanford H and 
ENS || experiments. With the current error bars, the 
ENS results do not show significant deviations from the 
classical 3D behavior. These results agree fairly well with 
our calculations. The Stanford experiment gives some- 
what lower values of Qth/^o at the lowest temperatures 
of the experiment. 
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FIG. 6. Thermalization rate versus temperature in the uni- 
tarity limit (a = 00). The solid curve shows the result of our 
calculations for Q t h/Qo, and the dotted line the 3D result 
Q,th /Oq. Squares and circles show the data of the Stanford 
and ENS experiments. 
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In the hydrodynamic regime for the gas cloud, where 
the characteristic collisional frequency greatly exceeds 
the radial frequency w, our assumption of quasiequilib- 
rium at instantaneous (time-dependent) values of T z and 
T p may not be valid. Nevertheless, the shape of the curve 
£l t h(T) qualitatively remains the same, including the ex- 
ponential decrease with temperature at T < Hujq and 
power law decrease with increasing T at temperatures 
larger than Tiloq. However, the maximum value of flth 
will be somewhat lower (in particular, of the order of lu 

§)• 

The number of particles (per "2D" sheet of atoms) in 
the Stanford experiment [|| was N ~ 10 4 p3| ], which is 
by a factor of 20 higher than at ENS for T m Tujq ||. We 
then estimate the 2D density of atoms for these tempera- 
tures to be n ~ 2.5 x 10 8 cm -2 at Stanford (uj w 90 Hz), 
and n ~ 0.5 x 10 8 cmT 2 at ENS (w w 180 Hz). For these 
densities, the ratio of the collisional frequency fl in Fig. 5 
to the radial frequency is fl/oj ~ 0.3 in the ENS exper- 
iment, and Q/u) ~ 3 in the experiment at Stanford. At 
temperatures T > hug the density n and the ratio Q/u> 
are smaller in both experiments. We thus see that the 
ENS experiment |8| was in the collisionless regime, al- 
though rather close to the hydrodynamic regime at tem- 
peratures T w TujJq. For these temperatures, the Stanford 
experiment Q has already achieved the hydrodynamic 
regime, and this can explain the discrepancy between our 
calculations and the Stanford results in Fig. 6. 

VII. INELASTIC 2-BODY PROCESSES 

Inelastic scattering of atoms is also influenced by the 
tight axial confinement of the particle motion. In this 
Section we will consider the inelastic 2-body processes, 
such as spin relaxation, in which the internal states of 
colliding atoms are changing, and the released internal- 
state energy of the atoms is transferred to their kinetic 
energy. Our goal is to establish a relation between the 
inelastic rates in 3D and those in the (tightly) axially 
confined geometry. The analysis given below relies on 
two important conditions widely met for the 2-body spin 
relaxation 

i) The energy release per collision greatly exceeds the gas 
temperature and the frequency of the axial confinement. 
Accordingly, the inelastic transitions occur at compara- 
tively short interparticle distances ~ Ri„ which are much 
smaller than the characteristic de Broglie wavelength of 
particles. 

ii) Thc inelastic transitions are caused by weak (spin- 
dipole, spin-orbit, etc.) interatomic interactions and can 
be treated with perturbation theory. 

To first order in perturbation theory the amplitude of 
inelastic scattering, defined in the same way as in the 
previous Sections, is given by a general expression |25j 

hn{e) = -2 / ^{v)U int {v)^ f {v)d\. (46) 



Here ipi(r) and ipf(r) are the true wavefunctions of the 
initial and final states of the relative motion of colliding 
atoms, and Ui n t{v) is the (weak) interatomic potential 
responsible for inelastic transitions. This potential is the 
same as in the 3D case. The function i/jf is also the 
same as in 3D, since the relative energy in the final state 
is much larger than Kloq. Thus, the only difference of 
the amplitude /,„ ( [46|) from the amplitude of inelastic 
scattering in the 3D case is related to the form of the 
wavefunction ipi. 

The characteristic interatomic distance Rt n at which 
the inelastic transitions occur, satisfies the inequality 
Rin <C A e (see item ii). Therefore, we are in the ultra- 
cold limit similar to that determined by Eq.(||) in the case 
of elastic scattering, and the conditions qRi„ -C 1 and 
Rin ar e satisfied. The former ensures a dominant 

contribution of the s-wave (of the initial wavefunction 
tpi) to the scattering amplitude /,- n (jiff). Due to the con- 
dition Ri n <C Iq, at distances r ~ Ri n the wavefunction 
ipi has a three-dimensional character: ipi(r) oc V^-dM) 
where ip3D{r) is the wavefunction of the 3D relative mo- 
tion at zero energy. For r ^> R e we have ipapOr) = 
(1 — a/r), and in order to be consistent with Eq.(n4h we 
should write 

iP i (r)= V (e)<p l ,(0)i> 3D (r), (47) 

where the coefficient 7?(e) is given by Eq. (|l7|) , and v is 
the quantum number of the initial state of the relative 
motion in the axial harmonic potential Vh{z). 

In the 3D case, the amplitude fo of inelastic scattering 
at zero initial energy is determined by Eq.(^6|) with ipi 
replaced by V>3D. Hence, Eq.([47|) directly gives a relation 
between the two scattering amplitudes: 

fin = v(e)<Pv(P)fo- (48) 

Due to the high relative kinetic energy of particles in the 
final state of the inelastic channel, the density of final 
states in this channel is independent of the axial con- 
finement. Therefore, relying on Eq.(f48|) the mean rate 
constant on n of inelastic collisions in the axially confined 
geometry and the corresponding collisional frequency can 
be represented in the form 

diin = (\rj{e)\ 2 tpl(0)6(e - hu> v))a ; Q in = a in n, (49) 

where ao is the 3D inelastic rate constant at zero energy. 

Note that in the ultra-cold limit the 3D inelastic rate 
constant is temperature independent and equal to ao if 
the scattering length \a\ < R e . For \a\ ^S> R e , the wave- 
function of the relative motion in the region of inter- 
atomic interaction takes the form ipi(r) — ?73d'03d('"), 
where 1 773 ^ | 2 = (1 + p 2 a 2 )~ 1 and p is the 3D relative 
momentum of colliding particles (see, e.g. [Pq| ). Hence, 
for the inelastic rate constant we have (ao(l + j^a 2 )" 1 ) . 
In the presence of axial confinement, averaging the fre- 
quency of inelastic collisions over the (quantum) axial 
density profile U3d(z), we obtain 
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a 



(l+p 2 a 2 ) 



L 3D 



dz. 



(50) 



The density profile n^D^z) accounts for the discrete struc- 
ture of quantum levels in the axial confining potential and 
for the quantum spatial distribution of particles. There- 
fore, Eq.(^o|) gives the ordinary 3D result only at temper- 
atures T 3> tiujo, where U3d(z) becomes the Boltzmann 
distribution ub{z). 

We first analyze the influence of axial confinement on 
Qin (|9|) for the case where \a\ <C At or, equivalently 
\a\ <C Iq and <jr|o| <C 1. In this case we may put rj = 1 
at any T, except for extremely low temperatures in the 
quasi2D regime. Then Eq.(f49|) gives 



((pl(0))a o n 



aon 



2nh 



tanh 



1/2 



hujQ 
~~f~ 



(51) 



One can easily check that Eq. ( pl| ) coincides with Eq. ( J50| ) 
in which p\a\ <C 1. The reason for this coincidence 
is that, similarly to the case of elastic scattering de- 
scribed by Eq.(p9[), for r\ = 1 the scattering ampli- 
tude fi n (^|) is independent of the relative energy e. 
Hence, the inelastic rate is influenced by the axial con- 
finement only through the axial distribution of parti- 
cles. However, this influence is significant, in contrast 
to the case of elastic scattering under the same condi- 
tions (see Eq.(|30|)). Qualitatively, for qr\a\ <C 1 we have 
ftin ~ «o n 3D- At temperatures T <§C hajg, a charac- 
teristic value of the 3D density is n^o ~ n/lo and we 
obtain fij n ~ ctQn/lo. For T S> hu>o, the 3D density 
Use ~ nfjnwo/T) 1 ' 2 and hence the frequency of inelastic 
collisions is f2j„ ~ {aon/l^^UQ/T) 1 / 2 . 

We now discuss the temperature dependence of the in- 
elastic rate for the case where \a\ > Zo , which in the ultra- 
cold limit (0) assumes that \a\ ^> R e . In the quasi2D 
regime and in the temperature interval e* < T < fkuo of 
the confinement-dominated 3D regime, the most impor- 
tant contribution to the inelastic rate in Eq.(^) comes 
from collisions with the axial quantum number v = 0. 
Then, using Eq. ( |l8| ) we express the parameter rj through 
the elastic amplitude /oo and obtain a relation between 
fli n and the mean frequency of elastic collisions f2(T): 



n in (T) = (|/ 00 (e)| 5 



(4™^ (0)) 2 



!l(T)0, (52) 



where (3 = (l/1287r 3 ) 1 / 2 (mZ a /?ia 2 ) is a dimensionless 
parameter independent of temperature. The tempera- 
ture dependence of Q is displayed in Fig.l and Fig. 5 and 
was discussed in sections [fyjand |v|. Note that the param- 
eter (3 is not equal to zero for \a\ — ► oo. In this case, since 
the amplitude /o is calculated with the wavefunction tp3D 
which behaves as a/r for r — ► oo, we have ao oc a 2 and 
j3 = const. 

At temperatures T 3> hujQ, using the same method 
as in Section ^ for the case of elastic scattering, from 
Eq.(fl9|) we recover the 3D result Slf® given by Eq.(Ec"' 



with U3d(z) — ub{z). In the limiting case, where 



T 



1/2 



(53) 



Comparing Eq.(p3[) with Eq.(p2[), we see that in the 
confinement-dominated 3D regime the deviation of the 
inelastic rate from the ordinary 3D behavior should be 
larger than that in the case of elastic scattering. 

As follows from Eq.(|52|) and Fig. 5, for \a\ > Iq the 
inelastic frequency Oi„ reaches its maximum at tem- 
peratures near the border between the quasi2D and 
confinement-dominated 3D regimes. The maximum 
value of Qi„ is close to 



a, 



I6h 



0. 



(54) 



From Eq.([49|) one finds that at any T the ratio fim/fim 
depends only on two parameters: T/Hloq and a/lo. In 
Fig. 7 we present our numerical results for Q.i n /VLi n as a 
function of T/hioo for a/lo equal to — 1, 1, and oo. As 
expected, the deviations from the 3D behavior are the 
largest in the unitarity limit. 
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FIG. 7. The quantity Qi n /Qi„ versus temperature. In (a) 
the parameter a/lo = — 1 (dashed curve), and a/lo — 1 (dot- 
ted curve). In (b) a/lo = oo (unitarity limit). The solid 
curves in (a) and (b) show the 3D result (E2I) . 
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The inelastic rate of spin relaxation in a tightly (axi- 
ally) confined gas of atomic cesium has been measured for 
the unitarity limit in the Stanford experiment 0. Due 
to a shallow radial confinement of the cloud in this ex- 
periment, the 2D density n ~ 1/T (see Eq. (pq) ) . Then, 
Eq.© gives the 3D inelastic frequency Q?*~~ 1/T 5 / 2 , 
whereas in the temperature interval £* < T < Tlloq of the 
confinement-dominated regime Eqs. (|3^) and ( |52] ) lead 
to fi ira ~ (1/T)(1 — exp (—huo/T)). In order to compare 
our calculations with the data of the Stanford experi- 
ment on spin relaxation, in Fig. 8 we display the ratio of 
Qin(T '/Jtioq) to £li n at T = Shujo which was the highest 
temperature in the experiment. The temperature depen- 
dence of the inelastic rate, following from the Stanford 
results, agrees fairly well with the calculations and shows 
significant deviations from the 3D behavior. It should be 
noted that, in contrast to thermalization rates, the in- 
elastic decay rate is not sensitive to whether the gas is in 
the collisionless or hydrodynamic regime [0]. 




1 2 3 4 



T/hG) Q 

FIG. 8. The ratio of Qi n (T/fkJo) to Qi n at T = 3huio versus 
temperature in the unitarity limit (a = oo). Fhe solid curve 
shows the result of our numerical calculations, and the dot- 
ted line the 3D limit. Circles show the data of the Stanford 
experiment. 

VIII. CONCLUDING REMARKS 

In conclusion, we have developed a theory which de- 
scribes the influence of a tight axial confinement of the 
particle motion on the processes of elastic and inelastic 
scattering. The most interesting case turns out to be 
the one in which the 3D scattering length a exceeds the 
extension of the wavefunction in the axial direction, Iq. 
In the ultra-cold limit defined by Eq.([|), the condition 
\a\ > Iq automatically requires large \a\ compared to the 
radius of interatomic interaction R e . Then we have a 
pronounced confinement-dominated 3D regime of scat- 
tering at temperatures on the order of Hcoq. Treating 
interatomic collisions as three-dimensional, the relative 
momentum of colliding atoms is related to the quantum 



character of the axial motion in the confining potential 
and becomes of the order of 1/Iq. As a result, the scat- 
tering rate can strongly deviate from the ordinary 3D be- 
havior. The axial extension of the wavefunction, achieved 
in the experiments at Stanford and ENS |5|-|8| , is Iq ~ 200 
A. The required value of the scattering length, \a\ > Iq 
and \a\ 3> R ei is characteristic for Cs atoms (R e ~ 100A ) 
and can also be achieved for other alkali atoms by using 
Feshbach resonances. 

In order to observe the 2D features of the particle mo- 
tion in the rates of interatomic collisions one has to reach 
the quasi2D regime of scattering, which requires much 
lower temperatures, at least by an order of magnitude 
smaller than Hujq. For ujq ~ 80kHz (Huq ~ 4 /iK) as in 
the Stanford ||[| and ENS @,|| experiments, these are 
temperatures below 400nK. As one can see from Fig. 5, 
the rate of elastic collisions is still rather large for these 
temperatures and, hence, one can think of achieving them 
by evaporative cooling. Moreover, for realistic radial fre- 
quencies u> ~ 100Hz there is a hope to achieve quantum 
degeneracy and observe a cross-over to the BEC regime. 
The cross-over temperature is T c sa N^-^huj (see (27]] and 
the discussion in |^]), and for N ~ 1000 particles in a 
quasi2D layer we find T c sa lOOnK. 

Another approach to reach BEC in the quasi2D regime 
will be to prepare initially a 3D trapped condensate and 
then adiabatically slowly turn on the tight axial confine- 
ment. Manipulating the obtained (quasi)2D condensate 
and inducing the appearance of thermal clouds with tem- 
peratures T < T c , one can observe interesting phase co- 
herence phenomena originating from the phase fluctua- 
tions of the condensate in quasi2D (see ||). 

Interestingly, at temperatures T ~ T c the collisional 
frequency Q can be on the order of the cross-over tem- 
perature if \a\ > lo- This follows directly from Fig.l 
and Eq.([2^) which give fl ~ irfin/m even at tempera- 
tures by two orders of magnitude smaller than Hloq. As 
the 2D density of thermal particles is n ~ Nmui 2 /T, we 
immediately obtain HQ(T C ) ~ N 1 / 2 ^ sa T c . This con- 
dition means that the trapped gas becomes a strongly 
interacting system. The mean free path of a particle, 
v/Q(T c ), is already on the order of its de Broglie wave- 
length h/yJmT c . At the same time, the system remains 
dilute, since the mean interparticle separation is still 
much larger than the radius of interatomic interaction 
R e . In this respect, the situation is similar to the 3D 
case with a large scattering length a ^> R e at densities 
where na 3 ~ 1. The investigation of the cross-over to the 
BEC regime in such strongly interacting quasi2D gases 
should bring in analogies with condensed matter systems 
or dense 2D gases. Well-known examples of dense 2D 
systems in which the Kostcrlitz-Thouless phase transi- 
tion ps[ has been found experimentally, are monolayers 
of liquid helium |^{| and the quasi2D gas of atomic hy- 
drogen on liquid helium surface ]30[ . 

On the other hand, for \a\ <C Iq the collisional fre- 
quency near the BEC cross-over, Sl(T c ) <g; T c /h, and the 
(quasi)2D gas remains weakly interacting. Then, the na- 
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ture of the cross-over is questionable (see ||). Generally 
speaking, one can have both the ordinary BEC cross- 
over like in an ideal trapped gas |27]] and the Kosterlitz- 
Thouless type |2^] of a cross-over. We thus see that ax- 
ially confined Bose gases in the quasi2D regime are re- 
markable systems where by tuning a or Iq one can modify 
the nature of the BEC cross-over. 
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